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1 INTRODUCTION: 



Let us revisit briefly the foundation of n-dimensional elasticity theory as it can be found today 
in any textbook, restricting our study to n = 2 for simplicity. If x = (x^, x^) is a point in the plane 
Ph ■ and ^ = (^^(x), f^(x)) is the displacement vector, lowering the indices by means of the Euclidean 

I metric, we may introduce the "small" deformation tensor e = (etj = eij = {l/2){di£_j + dj£,i)) with 

n{n + l)/2 — 3 (independent) components (eii,ei2 — £21, ^22)- If we study a part of a deformed 
body, for example a thin elastic plane sheet, by means of a variational principle, we may introduce 
the local density of free energy (p{e) — ip{eij\i < j) = i^(eii, £12, £22) and vary the total free energy 
F — Jip{e)dx with dx = dx^ A dx^ by introducing cr*^ — dip/deij for i < j in order to obtain 
SF = J{a^^5€ii + (T^^Sei2 + a-^^6€22)dx. Accordingly, the "decision" to define the stress tensor 
' (7 by a symmetric matrix with = is purely artificial within such a variational principle. 

. Indeed, the usual Cauchy device (1828) assumes that each element of a boundary surface is acted 

on by a surface density of force a with a linear dependence a = {a^^{x)nr) on the outward normal 
unit vector ft = (ur) and does not make any assumption on the stress tensor. It is only by an equi- 
librium of forces and couples, namely the well known phenomenological static torsor equilibrium,, 
that one can "prove" the symmetry of cr. However, even if we assume this symmetry, we now need 
' the different summation a^^ Seij = a^^Sen + 2a^^Sei2 + a^^6e22 = cr^^drS^i. An integration by part 

and a change of sign produce the volume integral Jdr{cr^^)S£,idx leading to the stress equations 
drcr^^ — 0. The classical approach to elasticity theory, based on invariant theory with respect to the 
group of rigid motions, cannot therefore describe equilibrium of torsors by means of a variational 
principle where the proper torsor concept is totally lacking. 



There is another equivalent procedure dealing with a variational calculus with constraint. In- 
deed, as we shall see in Section 7, the deformation tensor is not any symmetric tensor as it must 
satisfy n^(ri^ — 1)/12 compatibility conditions (CC), that is only 922£ii + 9ii£22 — 29i2£i2 = 
when n — 2. In this case, introducing the Lagrange multiplier —(f) for convenience, we have to 
vary J{(p{e) — (f>{d22^ii +9ii£22 — 29i2£i2))dx for an arbitrary e. A double integration by parts now 
provides the parametrization a^^ = 9220,0"^^ = a'^'^ = — 9i20, cr^^ = dii4> of the stress equations 
by means of the Airy function and the formal adjoint of the CC, on the condition to observe that 
we have in fact 2a^^ ~ — 29i20 as another way to understand the deep meaning of the factor "2" 
in the summation. In arbitrary dimension, it just remains to notice that the above compatibility 
conditions arc nothing else but the linearized Riemann tensor in Riemannan geometry, a crucial 
mathematical tool in the theory of general relativity. 

It follows that the only possibility to revisit the foundations of engineering and mathematical 
physics is to use new mathematical methods, namely the theory of systems of partial differential 
equations and Lie pseudogroups developped by D.C. Spencer and coworkers during the period 
1960-1975. In particular. Spencer invented the first order operator now wearing his name in order 
to bring in a canonical way the formal study of systems of ordinary differential (OD) or partial dif- 
ferential (PD) equations to that of equivalent first order systems. However, despite its importance. 
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the Spencer operator is rarely used in mathematics today and, up to our knowledge, has never 
been used in engineering or mathematical physics. The main reason for such a situation is that 
the existing papers, largely based on hand-written lecture notes given by Spencer to his students 
(the author was among them in 1969) are quite technical and the problem also lies in the only 
"accessible" book "Lie equations" he published in 1972 with A. Kumpera. Indeed, the reader can 
easily check by himself that the core of this hook has nothing to do with its introduction recalling 
known differential geometric concepts on which most of physics is based today. 

The first and technical purpose of this chapter, an extended version of a lecture at the second 
workshop on Differential Equations by Algebraic Methods (DEAM2, february 9-11, 2011, Linz, 
Austria), is to recall briefly its definition, both in the framework of systems of linear ordinary or 
partial differential equations and in the framework of differential modules. The local theory of Lie 
pseudogroups and the corresponding non-linear framework are also presented for the first time in 
a rather elementary manner though it is a difficult task. 

The second and central purpose is to prove that the use of the Spencer operator constitutes 
the common secret of the three following famous books published about at the same time in the 
beginning of the last century, though they do not seem to have anything in common at first sight 
as they are successively dealing with the foundations of elasticity theory, commutative algebra, 
electromagnetisni (EM) and general relativity (GR): 

[C] E. and F. COSSERAT: "Theorie des Corps Deformables" , Hermann, Paris, 1909. 
[M] F.S. MACAULAY: "The Algebraic Theory of Modular Systems", Cambridge, 1916. 
[W] H. WEYL: "Space, Time, Matter", Springer, Berlin, 1918 (1922, 1958; Dover, 1952). 

Meanwhile we shall point out the striking importance of the second book for studying iden- 
tifiability in control theory. We shall also obtain from the previous results the group theoretical 
unification of finite elements in engineering sciences (elasticity, heat, electromagnetism) , solving 
the torsor problem and recovering in a purely mathematical way known field-matter coupling phe- 
nomena (piezzoelectricity, photoelasticity, streaming birefringence, viscosity, ...). 

As a byproduct and though disturbing it may be, the third and perhaps essential purpose is 
to prove that these unavoidaMe new differential and homological methods contradict the existing 
mathematical foundations of both engineering (continuum mechanics, electromagnetism) and math- 
ematical (gauge theory, general relativity) physics. 

Many explicit examples will illustate this chapter which is deliberately written in a rather self- 
contained way to be accessible to a large audience, which does not mean that it is elementary in 
view of the number of new concepts that must be patched together. However, the reader must 
never forget that each formula appearing in this new general framework has been used explicitly 
or implicitly in [C], [M] and [W] for a mechanical, mathematical or physical purpose. 

2 FROM LIE GROUPS TO LIE PSEUDOGROUPS 

Evariste Galois (1811-1832) introduced the word "group" for the first time in 1830. Then the 
group concept slowly passed from algebra (groups of permutations) to geometry (groups of trans- 
formations). It is only in 1880 that Sophus Lie (1842-1899) studied the groups of transformations 
depending on a finite number of parameters and now called Lie groups of transform,ations. Let X 
be a manifold with local coordinates x = {x^, ...,x") and G be a Lie group, that is another manifold 
with local coordinates a = (a^, a^) called parameters with a composition GxG G : {a,b) ^ ab, 
an inverse G ^ G : a ^ and an identity e G G satisfying: 

{ab)c = a{bc) = abc, aa~^ = a~^a = e, ae = ea = a, Va, b,c G G 

DEFINITION 2.1: G is said to act on X if there is a map X xG ^ X : {x,a) ^ y = ax = f{x,a) 
such that {ah)x = a{bx) = abx,\/a,b £ G,Vx G X and, for simplifying the notations, we shall use 
global notations even if only local actions are existing. The set Gx = {a £ G \ ax = x} is called 
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the isotropy subgroup of G at x <E X. The action is said to be effective if ax = x, Vx S X a = e. 
A subset C X is said to be invariant under the action of G if aS C S*, Va G G and the orbit of 
X € X is the invariant subset Gx = {ax \ a ^ G} C X. If G acts on two manifolds X and Y, a 
map / : X — ^ y is said to be equivariant if f{ax) = af{x),Mx G X,Va e G. 

For reasons that will become clear later on, it is often convenient to introduce the graph 
X G ^ X X X : {x,a) {x,y = ax) of the action. In the product X x X, the first fac- 
tor is called the source while the second factor is called the target. 

DEFINITION 2.2: The action is said to be free if the graph is injective and transitive if the 
graph is surjective. The action is said to be simply transitive if the graph is an isomorphism and 
X is said to be a principal homogeneous space (PHS) for G. 

In order to fix the notations, we quote without any proof the " Three Fundamental Theorems 
of Lie" that will be of constant use in the sequel ([26]): 

FIRST FUNDAMENTAL THEOREM 2.3: The orbits x = f{xo,a) satisfy the system of 
PD equations dx^/da'^ = 9p{x)LuP{a) with det{u;) ^ 0. The vector fields 9p = 9p{x)dt are called 
infinitesimal generators of the action and are linearly independent over the constants when the 
action is effective. 

If X is a manifold, we denote as usual by T = T{X) the tangent bundle of X, by T* = T*{X) 
the cotangent bundle, by A'^T* the bundle of r-forms and by SqT* the bundle of q-symmetric ten- 
sors. More generally, let £ he a fibered manifold, that is a manifold with local coordinates (.t*, y'^) 
for i = l,...,n and k = l,...,m simply denoted by {x,y), projection n : £ ^ X : {x,y) ^ (x) and 
changes of local coordinates x = ifi{x),y = ^p{x,y). If £ and are two fibered manifolds over X 
with respective local coordinates (x.y) and {x,z), we denote by £xxJ' the fibered product of £ 
and J" over X as the new fibered manifold over X with local coordinates {x,y,z). We denote by 
/ : X — )• £ : (x) — > (.X, y = /(x)) a global section of £, that is a map such that tt o / = idx but 
local sections over an open set U C X may also be considered when needed. Under a change of 
coordinates, a section transforms like f{ip{x)) = tp{x, /(x)) and the derivatives transform like: 

^(^(xM^^ix) = ^{x,f{x)) + ^{x,f{xmf\x) 
We may introduce new coordinates {x'^ ,y^ ,yf) transforming like: 

yidiip^ix) = ■^ix,y) + ^{x,y)yf 

We shall denote by Jq{£) the q-jet bundle of £ with local coordinates {x^,y'',yf,yfj, ...) = {x,yq) 
called jef coordinates and sections fq : (x) (x, /'^(x), /f (x), /j^ (x), ...) = (x, /g(x)) transforming 
like the sections jq{f) : (x) {x,f^{x),dif^{x),dijf^{x),...) = {x, jq{f){x)) where both fq and 
jq{f) are over the section / of £. Of course Jq{£) is a fibered manifold over X with projection tt, 
while Jq+r{£) is a fibered manifold over Jq{£) with projection 7r^+'",Vr > 0. 

DEFINITION 2.4: A system of order g on £" is a fibered submanifold TZq C Jq{£) and a solution 
of TZq is a section f of £ such that jq{f) is a section of TZq. 

DEFINITION 2. 5: When the changes of coordinates have the linear form x = 1^(2;), y — A{x)y, 
we say that £ is a vector bundle over X and denote for simplicity a vector bundle and its 
set of sections by the same capital letter E. When the changes of coordinates have the form 
X = <p{x), y = A(x)y + B{x) we say that £ is an affine bundle over X and we define the associated 
vector bundle E over X by the local coordinates {x,v) changing like x = ^{x),v = A{x)v. 

DEFINITION 2.6: If the tangent bundle T{£) has local coordinates {x,y,u,v) changing like 
= dufP {x)u^ ,v'' = ^^(x, + |^(x, t/)u^, we may introduce the vertical bundle V{£) C T{£) 
as a vector bundle over £ with local coordinates (x, y, v) obtained by setting u = and changes 
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= fpr(3^! y)v''- Of course, when £ is an afEne bundle with associated vector bundle E over X, 
we have V{£) = £ xx E. 

For a later use, if f is a fibcrcd manifold over X and / is a section of £, wc denote by f~^{V{£)) 
the reciprocal image of V^(f) by / as the vector bundle over X obtained when replacing {x, y, v) 
by {x,f{x),v) in each chart. It is important to notice in variational calculus that a variation 6f 
of / is such that Sf{x), as a vertical vector field not necessary "small", is a section of this vector 
bundle and that (/, (5/) is nothing else than a section of V{£) over X. 

We now recall a few basic geometric concepts that will be constantly used. First of all, if ^, ?7 G 
T, we define their bracket 77] G T by the local formula ([^, ri]y{x) = £J'{x)drrf{x) — rj^{x)dsC{^) 
leading to the Jacobi identity [r?, C]] + [v^ [C)C]] + [Cj [Cj^]] = 0)V^)??)C G T allowing to define a 
Lie algebra and to the useful formula T(/)(r?)] = T(/)([^, 77]) where T(/) : T{X) T{Y) 

is the tangent mapping of a map f : X ^Y. 

SECOND FUNDAMENTAL THEOREM 2.7: If 61, Op are the infinitesimal generators of 

the effective action of a lie group G on X, then [9p,da-] = c^^Or where the c^^ are the structure 
constants of a Lie algebra of vector fields which can be identified with Q = TgiG). 

When / = {ii < ... < v} is a multi- index, we may set dx^ = dx*^ A ... A dx^"^ for describing 
A^T* and introduce the exterior derivative d : A^T* — > /\^+'^T* : u> = LVjdx^ doj = diLVjdx'^ Adx^ 
with d'^ = do d = in the Poincare sequence: 

^ ^ A"'"T'* ^ > A^r* ^ ) ^ V, ^ Q 

The Lie derivative of an r-form with respect to a vector field ^ € T is the linear first order 
operator jC{^) linearly depending on ji(^) and uniquely defined by the following three properties: 

1) mf = e/ = edif,yf G A^T* = c-(x). 

2) Ci^d = dC{i). 

3) C{0{a A 13) = mOa) A /3 + a A (£(0/9), Va, /3 G AT*. 

It can be proved that £(0 = i(0(i + di{(_) where i{(_) is the interior multiplication {i{^)u!)ii = 
eoJu,...i. and that [£(?), £(r/)] = m o £{ri) - C{v) o £{^) = £([e, r?]), V^, 7? G T. 

Using crossed-derivatives in the PD equations of the First Fundamental Theorem and introduc- 
ing the family of 1-forms w"^ = u}'^{a)da'^ both with the matrix a = of right invariant vector 
fields, we obtain the compatibility conditions (CC) expressed by the following corollary where one 
must care about the sign used: 

COROLLARY 2.8: One has the Maurer-Cartan (MC) equations dio'' + cJ,^u}P Alj" = or the 
equivalent relations [Qp,aCT] = cj,^ar- 

Applying d to the MC equations and substituting, we obtain the integrability conditions (IC): 

THIRD FUNDAMENTAL THEOREM 2.9: For any Lie algebra Q defined by structure con- 
stants satisfying : 

c'' ^0 + c^ c'' +c^ c^ =0 

one; can construct an analytic group G such that Q = Te{G). 

EXAMPLE 2.10: Considering the affine group of transformations of the real line y = a^x + a^, 
we obtain 9i — xdx,02 = ^ [^1)^2] = —O2 and thus uj^ = {l/a^)da^ ,uj'^ — da^ — {a'^/a^)da^ ^ 
dio^ = 0, rfw^ - A = 0. 

Only ten years later Lie discovered that the Lie groups of transformations are only particular 
examples of a wider class of groups of transformations along the following definition where aut{X) 
denotes the group of all local diffeomorphisms of X: 
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DEFINITION 2.11: A Lie pseudogroup of transformations F C aut{X) is a group of transfor- 
mations solutions of a system of OD or PD equations such that, if y = f{x) and z = g{y) are 
two solutions, called finite transformations, that can be composed, then z = g o f{x) = h{x) and 
X = f~^{y) = g{y) are also solutions while y = xis & solution. 

The underlying system may be nonlinear and of high order as we shall sec later on. We shall 
speak of an algebraic pseudogroup when the system is defined by differential polynomials that is 
polynomials in the derivatives. In the case of Lie groups of transformations the system is obtained 
by differentiating the action law y = f{x, a) with respect to x as many times as necessary in order 
to eliminate the parameters. Looking for transformations "close" to the identity, that is setting 
y = X + t^{x) + ... when t ^ 1 is a small constant parameter and passing to the limit f — >■ 0, we 
may linearize the above nonlinear system of finite Lie equations in order to obtain a linear system 
of infinitesimal Lie equations of the same order for vector fields. Such a system has the property 
that, if ^, T] are two solutions, then rj] is also a solution. Accordingly, the set C T of solutions 
of this new system satifies [O, G] C O and can therefore be considered as the Lie algebra of T. 

Though the collected works of Lie have been published by his student F. Engel at the end of 
the 19*'' century, these ideas did not attract a large audience because the fashion in Europe was 
analysis. Accordingly, at the beginning of the 20*'' century and for more than fifty years, only two 
frenchmen tried to go further in the direction pioneered by Lie, namely Elie Cartan (1869-1951) 
who is quite famous today and Ernest Vessiot (1865-1952) who is almost ignored today, each one 
deliberately ignoring the other during his life for a precise reason that we now explain with details 
as it will surprisingly constitute the heart of this chapter. (The author is indebted to Prof. Maurice 
Janet (1888-1983), who was a personal friend of Vessiot, for the many documents he gave him, 
partly published in [25]). Roughly, the idea of many people at that time was to extend the work 
of Galois along the following scheme of increasing difficulty: 

1) Galois theory ([34]): Algebraic equations and permutation groups. 

2) Picard- Vessiot theory ([17]): OD equations and Lie groups. 

3) Differential Galois theory ([24], [37]): PD equations and Lie pseudogroups. 

In 1898 Jules Drach (1871-1941) got and published a thesis ([9]) with a jury made by Gaston 
Darboux, Emile Picard and Henri Poincare, the best leading mathematicians of that time. How- 
ever, despite the fact that it contains ideas quite in advance on his time such as the concept of 
a "differential field" only introduced by J.-F. Ritt in 1930 ([31]), the jury did not notice that the 
main central result was wrong: Cartan found the counterexamples, Vessiot recognized the mistake, 
Paul Painleve told it to Picard who agreed but Drach never wanted to acknowledge this fact and 
was supported by the influent Emile Borel. As a byproduct, everybody flew out of this "affair", 
never touching again the Galois theory. After publishing a prize- winning paper in 1904 where he 
discovered for the first time that the differential Galois theory must be a theory of (irreducible) 
PHS for algebraic pseudogroups, Vessiot remained alone, trying during thirty years to correct the 
mistake of Drach that we finally corrected in 1983 ([24]). 

3 CARTAN VERSUS VESSIOT : THE STRUCTURE EQUATIONS 

We study first the work of Cartan which can be divided into two parts. The first part, for 
which he invented exterior calculus, may be considered as a tentative to extend the MC equations 
from Lie groups to Lie pseudogroups. The idea for that is to consider the system of order q and its 
prolongations obtained by differentiating the equations as a way to know certain derivatives called 
principal from all the other arbitrary ones called parametric in the sense of Janet ([13]). Replacing 
the derivatives by jet coordinates, we may try to copy the procedure leading to the MC equations 
by using a kind of "composition" and "inverse" on the jet coordinates. For example, coming back 
to the last definition, we get successively: 



Now if g = f ^ then g o f = id and thus = 1, ••. while the new identity idg = jq{id) is made 

by the successive derivatives of y = x, namely (1,0,0,...). These awfully complicated computa- 



dh dg df 
dx dy dx 
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tions bring a lot of structure constants and have been so much superseded by the work of Donald 
C. Spencer (1912-2001) ([11], [12], [18], [33]) that, in our opinion based on thirty years of explicit 
computations, this tentative has only been used for classification problems and is not useful for 
applications compared to the results of the next sections. In a single concluding sentence. Car- 
tan has not been able to "go down" to the base manifold X while Spencer did succeed fifty years later. 

We shall now describe the second part with more details as it has been (and still is !) the cru- 
cial tool used in both engineering (analytical and continuum mechanics) and mathematical (gauge 
theory and general relativity) physics in an absolutely contradictory manner. We shall try to use 
the least amount of mathematics in order to prepare the reader for the results presented in the 
next sections. For this let us start with an elementary experiment that will link at once continuum 
mechanics and gauge theory in an unusual way. Let us put a thin elastic rectilinear rubber band 
along the x axis on a flat table and translate it along itself. The band will remain identical as 
no deformation can be produced by this constant translation. However, if we move each point 
continuously along the same direction but in a point depending way, for example fixing one end 
and pulling on the other end, there will be of course a deformation of the elastic band according 
to the Hooke law. It remains to notice that a constant translation can be written in the form 
y = X + as in Example 2.10 while a point depending translation can be written in the form 
y = X + a'^{x) which is written in any textbook of continuum mechanics in the form y = x + ^{x) 
by introducing the displacement vector ^. However nobody could even imagine to make also 
point depending and to consider y = a^{x)x + a^(x) as we shall do in Example 7. 9. We also notice 
that the movement of a rigid body in space may be written in the form y = a{t)x + b{t) where now 
a{t) is a time depending orthogonal matrix and b{t) is a time depending vector. What makes all 
the difference between the two examples is that the group is acting on x in the first but not acting 
on t in the second. Finally, a point depending rotation or dilatation is not accessible to intuition 
and the general theory must be done in the following manner. 

If X is a manifold and G is a lie group not acting necessarily on X, let us consider maps 
a : X G : (x) {a{x)) or cquivalently sections of the trivial (principal) bundle XxG over X. 
If x + dx is a point of X close to x, then T(a) will provide a point a + da = a+^dx close to a on G. 
We may bring a back to e on G by acting on a with o~^, either on the left or on the right, getting 
therefore a 1-form a~^da = ^ or daa~^ = B. As aa~^ = e we also get daa~^ = —ada~^ = —b~^db 
if we set b = a~^ as a way to link A with B. When there is an action y — ax, we have x = a~^y = by 
and thus dy = dax = daa~^y, a result leading through the First Fundamental Theorem of Lie to 
the equivalent formulas: 

a-^da = A^ {Alix)dx' = -ojl{b{x))d^b'' {x)dx') 

daa-^ = B^ {BJ{x)dx^ = ujUa{x))d^a'' {x)dx') 

Introducing the induced bracket [^,A](^,ry) = [74(^), A(ry)] e g,y^,r] £ T, we may define the 2- 
form dA - [A, A] = F G A^T* (g) by the local formula (care to the sign): 

d,A^x) - d,Al{x) ~ c;^A^{x)A^x) = Flj{x) 
and obtain from the second fundamental theorem: 

THEOREM 3.1: There is a nonlinear gauge sequence: 

XxG ^ T*®g ^ A^T*®^ 
a — > a-'^da = A — > dA-[A,A]=F 

Choosing a "close" to e, that is a{x) = e + tX{x) + ... and linearizing as usual, we obtain the 
linear operator d : A°T* ®g ^ A^T* ® Q : {X'ix)) (5iA'"(a;)) leading to: 

COROLLARY 3.2: There is a linear gauge sequence: 

A°T* (g) ^? A A^T* ®g-^ A^T* ® a A ... A A"T* ® a — ^ 
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which is the tensor product by Q of the Poincare sequence: 

REMARK 3.3: When the physicists C.N. Yang and R.L. Mills created (non-abelian) gauge the- 
ory in 1954 ([38], [39]), their work was based on these results which were the only ones known at 
that time, the best mathematical reference being the well known book by S. Kobayashi and K. 
Nomizu ([15]). It follows that the only possihility to describe elccromagnetism (EM) within this 
framework was to call A the Yang-Mills potential and F the Yang-Mills field (a reason for choosing 
such notations) on the condition to have dim{G) = 1 in the abelian situation c = and to use a 
Lagrangian depending on F = dA — [A, A] in the general case. Accordingly the idea was to select 
the unitary group U{1), namely the unit circle in the complex plane with Lie algebra the tangent 
line to this circle at the unity (1,0). It is however important to notice that the resulting Maxwell 
equations dF = have no equivalent in the non-abelian case c^O. 

Just before Albert Einstein visited Paris in 1922, Cartan published many short Notes ([5]) 
announcing long papers ([6]) where he selected G to be the Lie group involved in the Poincare 
(conformal) group of space-time preserving (up to a function factor) the Minkowski metric U! = 
{dx^)"^ -\- {dx^Y + {dx^Y ~ {dx^Y with = ct where c is the speed of light. In the first case F is 
decomposed into two parts, the torsion as a 2-form with value in translations on one side and the 
curvature as a 2-form with value in rotations on the other side. This result was looking coherent at 
first sight with the Hilbert variational scheme of general relativity (GR) introduced by Einstein in 
1915 ([21], [38]) and leading to a Lagrangian depending on F = dA — [A, A\ as in the last remark. 

In the meantime, Poincare developped an invariant variational calculus ([22]) which has been 
used again without any quotation, successively by G. BirkhofF and V. Arnold (compare [4], 205-216 
with [2], 326, Th 2.1). A particular case is well known by any student in the analytical mechanics 
of rigid bodies. Indeed, using standard notations, the movement of a rigid body is described in 
a fixed Cartesian frame by the formula x{t) = a{t)xQ + b(t) where a{t) is a 3 x 3 time dependent 
orthogonal matrix (rotation) and b{t) a time depending vector (translation) as we already said. 
Differentiating with respect to time by using a dot, the absolute speed is v = x{t) = d{t)xo + b{t) 
and we obtain the relative speed a~^{t)v = {t)a{t)xi,) + ar^{t)b{t) by projection in a frame fixed 
in the body. Having in mind Example 2.10, it must be noticed that the so-called Eulerian speed 
V = v{x,t) = da~-^a;-|-6 — da~-^6 only depends on the 1-form B = {aa~^ ,b — da~^b). The Lagrangian 
(kinetic energy in this case) is thus a quadratic function of the 1-form A = {a~^d, a~^b) where 
a~^d is a 3 X 3 skew symmetric time depending matrix. Hence, "surprisingly'', this result is not 
coherent at all with EM where the Lagrangian is the quadratic expression {e/2)E'^ — (l/2/i)B^ 
because the electric field E and the magnetic field B are combined in the EM field as a 2-form 
satisfying the first set of Maxwell equations dF = 0. The dielectric constant e and the magnetic 
constant are leading to the electric induction D = eE and the magnetic induction H = (l//x)i? 
in the second set of Maxwell equations. In view of the existence of well known field-matter cou- 
plings such as piezoelectricity and photoelasticity that will be described later on, such a situation 
is contradictory as it should lead to put on equal footing 1-forms and 2-forms contrary to any 
unifying mathematical scheme but no other substitute could have been provided at that time. 

Let us now turn to the other way proposed by Vessiot in 1903 ([36]) and 1904 ([37]). Our 
purpose is only to sketch the main results that we have obtained in many books ([23-26], wc; do not 
know other references) and to illustrate them by a series of specific examples, asking the reader to 
imagine any link with what has been said. 

1) If £" = X X X, we shall denote by H, = Iiq{X,X) the open subfibered manifold of 
Jq{X X X) defined independently of the coordinate system by det{y^) ^ with source projec- 
tion aq : liq ^ X : {x,yq) — > (x) and target projection fiq : Wq ^ X : {x,yq) (y). We shall 
sometimes introduce a copy y of X with local coordinates {y) in order to avoid any confusion 
between the source and the target manifolds. Let us start with a Lie pseudogroup F c aut{X) 
defined by a system TZq C Tlq of order q. In all the sequel we shall suppose that the system is 
involutive (see next section) and that F is transitive that is Vx,y £ X,3f £T,y = f{x) or, equiv- 
alently, the map (ag,/?^) : TZq ^ X x X : {x,yq) ^ {x,y) is surjective. 
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2) The lie algebra C T of infinitesimal transformations is then obtained by linearization, 
setting y = X + t£,{x) + ... and passing to the limit t ^ in order to obtain the linear involutive 
system Rg = id~^ {V {TZg)) C Jq{T) by reciprocal image with 6 = e T\jg{$,) e Rg}. 

3) Passing from source to target, we may prolong the vortical infinitesimal transformations 
r] = 7]'' {y) to the jet coordinates up to order q in order to obtain: 

' ^'^'dy'^ dy-^' dyf ^dy^dy'"'^' dy-'^'^'dyf^ 

where we have replaced jg{f){x) by j/g, each component beeing the "formal" derivative of the 
previous one . 

4) As [0, 8] C 0, we may use the Frobenius theorem in order to find a generating fundamental 
set of differential invariants {^'^{yq)} up to order q which are such that ^^{yq) ~ ^^{Uq) by using 
the chain rule for derivatives whenever y = g{y) S F acting now on Y . Of course, in actual practice 
one must use sections of Rq instead of solutions but it is only in section 6 that we shall see why 
the use of the Spencer operator will be crucial for this purpose. Specializing the at idq {x) we 
obtain the Lie form ^'^{Vg) = u)^{x) of TZg. 

5) The main discovery of Vcssiot, fifty years in advance, has been to notice that the prolongation 
at order q of any horizontal vector field C = i^) commutes with the prolongation at order q 
of any vertical vector field r] = Vi^{y)-^, exchanging therefore the differential invariants. Keeping 
in mind the well known property of the Jacobian determinant while passing to the finite point of 
view, any (local) transformation y = f{x) can be lifted to a (local) transformation of the differen- 
tial invariants between themselves of the form u X{u,jg{f){x)) allowing to introduce a natural 
bundle T over X by patching changes of coordinates x = (p{x),u = X{u,jq{(p){x)). A section cj of 

is called a geometric object or structure on X and transforms like uj{f{x)) — X{uj{x),jg{f){x)) 
or simply u) = jg{f){ijo). This is a way to generalize vectors and tensors {q = 1) or even con- 
nections {q = 2). As a byproduct we have F = {/ G aut{X)\^^{jq{f j) = jq[f)^^{uj) = w} 
as a new way to write out the Lie form and we may say that F preserves lo. We also obtain 
T^q = {fq G ^q\fq^{^) = ^} ■ Coming back to the infinitesimal point of view and setting 
ft = exp{t£) G aut{X),WS, G T, we may define the ordinary Lie derivative with value in LJ~-^{V{Tj) 
by the formula : 

= 2?„C = = j^Mftr\^)\t=o ^Q^iCe T\C{Ou; = 0} 

while we have x ^ x + t^{x) + ... w"^ — >■ m"^ + td^^'^ Lj.^ {u) + ... where fj, — (/xi, ...,iJLn) is a 
multi-index as a way to write down the system of infinitesimal Lie equations in the Medolaghi form: 

il^ = mOioy = -Ll^'{u:{x))d^e + Cdrio^ix) = 

7) By analogy with "special" and "general" relativity, we shall call the given section special 
and any other arbitrary section general. The problem is now to study the formal properties of the 
linear system just obtained with coefficients only depending on ji{uj), exactly like L.P. Eisenhart 
did for J' ~ S2T* when finding the constant Riemann curvature condition for a metric ui with 
det{u}) ^ ([26], Example 10, p 249). Indeed, if any expression involving uj and its derivatives is a 
scalar object, it must reduce to a constant because F is assumed to be transitive and thus cannot 
be defined by any zero order equation. Now one can prove that the CC for thus for uj too, 
only depend on the $ and take the quasi-linear symbolic form v = I{ui) = A{u)ux + B{u) ~ 0, 
allowing to define an affine subfibered manifold Bi C Ji [T) over F. Now, if one has two sections 
w and Q oi the equivalence problem is to look for / e aut{X) such that jq{f)^^(oj) = Q. When 
the two sections satisfy the same CC, the problem is sometimes locally possible (Lie groups of 
transformations, Darboux problem in analytical mechanics,...) but sometimes not ([23], p 333). 

8) Instead of the CC for the equivalence problem, let us look for the integrability conditions 
(IC) for the system of infinitesimal Lie equations and suppose that, for the given section, all the 
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cqTiations of order q + r are obtained by differentiating r times only the equations of order q, then 
it was claimed by Vessiot ([36] with no proof, see [26], p 209) that such a property is held if and 
only if there is an equivariant section c : — t- J*! : (x, w) — >■ (x, u,v = c{u)) where Ti = Ji{T)/B\ 
is a natural vector bundle over F with local coordinates {x,u,v). Moreover, any such equivariant 
section depends on a finite number of constants c called structure constants and the IC for the 
Vessiot structure equations I{u\) = c{u) are of a polynomial form J(c) = 0. 

9) Finally, when Y is no longer a copy of X, a system Aq C Jq{X x Y) is said to be an au- 
tomorphic system for a Lie pseudogroup T C aut{Y) if, whenever y = f{x) and y = f{x) are two 
solutions, then there exists one and only one transformation y = g{y) G F such that / = g o f. 
Explicit tests for checking such a property formally have been given in [24] and can be implemented 
on computer in the differential algebraic framework. 

EXAMPLE 3.4: {Principal homogeneous structure) When F is made by the translations = 
a;' + a*, the Lie form is $*^(t/i) = y^ = (Kronecker symbol) and the linearization is 9^^*^ = 0. 
The natural bundle is = T*Xx---XxT* (n times) with det{u}) ^ and the general Medolaghi 
form is oj^di^^ + £J'dr(jJl = <S=> [^,aT] = with r = 1, if a = (a^) = (jO~^. Using crossed 
drivatives, one finally gets the zero order equations: 

CdMp{x)ciix){d,Lo]ix) d,u:l{x))) = 
leading therefore (up to sign) to the n'^{n — l)/2 Vessiot structure equations: 

diUj]{x) - djulix) ^ c;;^wf (x)cjj(a;) 

This result proves that the MC equations are only examples of the Vessiot structure equations. We 
finally explain the name given to this structure ([26], p 296). Indeed, when X is a PHS for a lie group 
G, the graph of the action is an isomorphism and we obtain a map X x X ^ G : {x, y) — ^ (a(.x, y)) 
leading to a first order system of finite Lie equations yx = gi^ {x , a{x , y)) . In order to produce a 
Lie form, let us first notice that the general solution of the system of infinitesimal equations is 
^ = X^Ot with A = est. Introducing the inverse matrix (cj) = (luI) of the reciprocal distribution 
a = {ar} made by n vectors commuting with {9r}, we obtain X = est [^, a] = <^ ^{0'^ = 0- 

EXAMPLE 3.5: [Ajfine and projective structures of the real line) In Example 2.10 with n = 1, 
the special Lie equations are ^■(^2) = yxx/Vx = dxx^ = with q = 2 and we let the reader 
check as an exercise that the general Lie equations are: 

— + ^{y)yx = oj{x) ^ dxxi + '^{x)dxi + idxuj{x) = 

Vx 

with no IC. The special section \s uj{x) =0. 

We could study in the same way the group of projective transformations of the real line y = 
{ax + b)/ {cx + d) and get with more work the general lie equations: 

— - li — f + ^{y)yl = ^{x) ^ dxxx^ + 2uj{x)dx^ + ^dxuj{x) = 
yx 2 yx 

There is an isomorphism Ji{J^aff) — J^aff^x^proj '■ ji{^) {'^^7 = ^xOJ — (l/2)w^). 

EXAMPLE 3.6: n = 2,q = l,V = {y^ = f{x^),y'^ = x'^/{df{x^)/dx^)} where / is an arbitrary 
invertible map. The involutive Lie form is: 

<f\y,) ^ yV, = x^ $^(yr) ^ yV^ = 0, <^'{y,) ^ ^ yly', - ylyj = 1 

We obtain = T*Xx/\^T* and w = {a, (3) where a is a 1-form and /3 is a 2-form with special 
section cj = {x^dx^,dx^ A dx^). It follows that da/ (3 is a well defined scalar because P ^ 0. The 
Vessiot structure equation is da = c/3 with a single structure constant c which cannot have any- 
thing to do with a Lie algebra. Considering the other section u) = {dx^,dx^ A dx^), we get c = 0. 
As c 7^ c, the equivalence problem ji(/)~^(w) = u) cannot even be solved formally. 
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EXAMPLE 3.7: (Symplectic structure) With n = 2p, g = 1 and J" = A^T*, let w be a closed 
2-form of maximum rank, that is cLu = 0, det{uj) ^ 0. The equivalence problem is nothing else than 
the; Darboux problem in analytical mechanics giving the possibility to write locally u = ^dp /\dq 
by using canonical conjugate coordinates {q,p) = {position, momentum). 

EXAMPLE 3.8: {Contact structure) With n = 3, q = 1, w = dx^ - x^dx^ =^ w A dw = 
dx^ A dx"^ A dx^, let us consider F = {/ e aut{X)\ji{f)~^ {w) = pw}. This is not a Lie form 
but we get: 

ji{f)-'^{dw) =dji{f)-^{w) = pdw + dpAw^ ji{f)-^{w Adw) ^ p^{w Adw) 

The corresponding geometric object is thus made by a 1-form density u = {(Ji,lJ2,ujs) that trans- 
forms like a 1-form up to the division by the square root of the Jacobian determinant. The unusual 
general Medolaghi form is: 

= COr{x)d,C - {l/2)cO^{x)drC + C drCOi{x) = 

In a symbolic way w A is now a scalar and the only Vessiot structure equation is: 

U)l{d2UJ3 ~ + UJ2{d3UJl - diCJs) + W3(9iW2 — d2Uli) = c 

For the special section u = (1, —x^,0) we have c = 1. If we choose ui = (1,0,0) we may define 
r by the system y| = O^y^ — 0,y|y| — — yl but now c = and the equivalence problem 
ji{f)~^{uj) = u) cannot even be solved formally. These results can be extended to an arbitrary odd 
dimension with much more work ([24], p 684). 

EXAMPLE 3.9: {Screw and complex structures) {n — 2, q = 1) In 1878 Clifford introduced ab- 
stract numbers of the form x^ + e.x^ with = in order to study helicoidal movements in the 
mechanics of rigid bodies. We may try to define functions of these numbers for which a derivative 
may have a meaning. Thus, if f{x^ + ex^) = f^{x^,x'^) + ep{x^ ,x'^), then we should get: 

df ^{A + eB){dx^ + edx^) = Adx^ + e{Bdx^ + Adx'^) = df + edf 

Accordingly, we have to look for transformations = f^{x^,x'^),y'^ = p{x^,x^) satisfying the 
first order involutive system of finite Lie equations 2/2=0) 2/1 ~ 2/i = with no CC. As we have 
an algebraic Lie pseudogroup, a tricky computation ([24], p 467) allows to prove that F is made 
by the transformations preserving a mixed tensor with square equal to zero as follows: 

v\ v\y^ ( ^ \ / y\ y\\^(Q 

We get the Lie form = yl/yl = 0, $^ = {yl)^/ {y\y2 — ulvi) = 1 and let the reader exhibit F. 
Finally, introducing similarly the abstract number i such that = —1, we get the Cauchy-Riemann 
system yl ~ 2/i = 0; 2/2 + 2/i = with no CC defining complex analytic transformations and the 
correponding geometric object or complex structure is a mixed tensor with square equal to minus 
the 2x2 identity matrix as we have now: 

yi yiy'fo -^\fyl _ ^ o -1 
yl yi ) [1 J[yf yl 

EXAMPLE 3.10: {Riemann structure) If w is a section of 7" = S2T* with det{uj) 7^ we get: 
Lie form ^ij{yi) = ^ki{y)y'^y\ = ^ij{x) 

Medolaghi form = {C{£_)uj)ij = ojrj{x)di(J' + bJir{x)dj^'^ + ^'^ drOJij{x) = 

also called Killing system, for historical reasons. A special section could be the Euclidean metric 
when n= 1, 2, 3 as in elasticity theory or the Minkowski metric when n = 4 as in special relativity. 
The main problem is that this system is not involutive unless we prolong the system to order two 
by differentiating once the equations. For such a purpose, introducing uj~^ = (w*^) as usual, we 
may define: 

Christoffel symbols lij{x) = ^i^''^ {x){diiJrj{x) + djOJri{x) — drCOij{x)) = 7ji(x) 

This is a new geometric object of order 2 allowing to obtain, as in Example 3.5, an isomorphism 
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ji{u}) ~ (w,7) and the second order equations with ^ = gi: 
Lie form {y\j + 7^^ {yjyly -) = 7y- (a;) 

Medolaghi form . = (£(07)f,- = dij^ + l^rj{x)di^'' + liM^j^'' ' Yij{x)dre + Cdrl^x) = 

where (r*j ) is a section of S'2T* 05 T. Surprisingly, the following expression: 

Riemann tensor pf^j {x) = d^f^ {x) - {x) + 7[^. (a;)7^i (x) - 7J; (x)7^^- (a;) 

is still a first order geometric object and even a tensor as a section of A^T* ®T* ®T satisfying the 

purely algebraic relations : 

Accordingly, the IC must express that the new first order equations {^{Op)Hj = only linear 
combinations of the previous ones and we get the Vessiot structure equations: 

Puji^) = c{S^^ij{x) - S'^uuix)) 
describing the constant Riemannian curvature condition of Eisenhart [10]. Finally, as we have 
Prij{x) = di'ylj{x) — dj'yl^{x) = 0, we can only introduce the Ricci tensor pij{x) = plrj{x) = Pji{x) 
by contracting indices and the scalar curvature p{x) = lo^^ (x) pij (x) in order to obtain p(x) = 
n{n — l)c. It remains to obtain all these results in a purely formal way, for example to prove that 
the number of components of the Riemann tensor is equal to n^(n^ — 1)/12 without dealing with 
indices. 

REMARK 3.11: Comparing the various Vessiot structure equations containing structure con- 
stants, we discover at once that the many c appearing in the MC equations are absolutely on equal 
footing with the only c appearing in the other examples. As their factors are either constant, linear 
or quadratic, any identification of the quadratic terms appearing in the Riemann tensor with the 
quadratic terms appearing in the MC equations is definitively not correct or, in an equivalent but 
more abrupt way, the Cartan structure equations have nothing to do with the Vessiot structure 
equations. As we shall see, most of mathematical physics today is based on such a confusion. 

REMARK 3.12: Let us consider again Example 3.5 with dxxf{x)/dxf{x) = cj(x) and intro- 
duce a variation r]{f{x)) = Sf{x) as in analytical or continuum mechanics. We get similarly 
Sdxf = dxSf = ^dxf and so on, a result leading to Suj{x) = dxfC{r])uj{f{x)) where the Lie 
derivative involved is computed over the target. Let us now pass from the target to the source by 
introducing r] = £,dxf f^^x/ = dx^dxf + £,dxxf and so on, a result leading to the particu- 
larly simple variation 5il) = C{$_)Q over the soure. As another example of this general variational 
procedure, let us compare with the similar variations on which classical finite elasticity theory is 
based. Starting now with 0Jki{fix))dif'^{x)djf''(x) = u)ij{x), where w is the Euclidean metric, we 
obtain {Sui)ij{x) = dif^{x)djf^{x){C{rj)oj)ki{f{x?}) where the Lie derivative involved is computed 
over the target. Passing now from the target to the source as before, we find the particularly simple 
variation 5u) = C{^)u> over the source. For "small" deformations, source and target are of course 
identified but it is not true that the infinitesimal deformation tensor is in general the limit of the 
finite deformation tensor (for a counterexample, see [25], p 70). 

Introducing a copy K of X in the general framework, (/, Sf) must be considered as a section 
of V{X xY) = {X X Y)xyT{Y) = X X T{Y) over X. When / is invertible (care), then we 
may consider the map f : X ^ Y : (x) ^ {y = f{x)) and define ^ € T by ry = T{f){^) or 
rather 77 = ji(/)(C) in the language of geometric object, as a way to identify f~^{V{X x Y)) with 
T = T{X). When / = id, this identification is canonical by considering vertical vectors along the 
diagonal A = {{x, y) £ X x Y\y = x} and we get Sco = SI £ Fq = uj~^{V{T)). We point out that 
the above vertical procedure is a nice tool for studying nonlinear systems ([26], III, C and [27], III, 2). 

4 JANET VERSUS SPENCER : THE LINEAR SEQUENCES 

Let p = {pi, pn) be a multi-index with length \p\ = pi + ... + pn, class i'li pi = ... = pi-i = 
0,Pij^O and p + li = {pi, ...,pi_i,pi + l,pi+i, ...,Pn)- We set yq = {?/^]l <k <m,Q <\p\ < q) 
with j/^ = y^ when \p\ =0. If i? is a vector bundle over X with local coordinates {x'^,y^) for 
i = l,...,n and k = l,...,m, we denote by Jg{E) the q-jet bundle of E with local coordinates 
simply denoted by {x,yq) and sections fq : {x) — >■ {x, f'^{x), fi{x), fij{x), ...) transforming like 
the section jq{f) : (a;) {x,f''{x),dif''{x),dijf''{x),...) when / is an arbitrary section of E. 



11 



Then both fq G Jq{E) and jq{f) G Jq{E) arc over f £ E and the Spencer operator inst allows 
to distinguish them by introducing a kind of "difference" through the operator D : Jg+i(iJ) — > 
T* ® Jg{E) : /,+! ^ - fg+i with local components (dif^ix) - ff{x),dif^{x) - ftj{x),...) 

and more generally {D .fq+i)'fi = dif^{x) — Z^+i. (a;)- In a symbolic way, when changes of 
coordinates are not involved, it is sometimes useful to write down the components of D in the 
form di = di — Si and the restriction of D to the kernel Sg+iT* ^ E of the canonical projection 
7^1"^^ : Jq+i{E) Jq{E) is minus the Spencer map 5 = dx' A S^ : Sq+iT* ® E ^ T* ^ SqT* ® E. 
The kernel of D is made by sections such that fg+i = ji{fq) = j2{fq-i) = ■■■ = jq+i{f)- Fi- 
nally, if Rq C Jq{E) is a system of order q on E locally defined by linear equations ^'^{x,yq) = 
a1,^{x)y^ = and local coordinates (a;, z) for the parametric jets up to order g, the r- prolongation 
Rq+r = Pr{Rq) — Jr{Rq) H Jq+r{E) C Jr{Jq{E)) is locally defined when r = 1 by the lin- 
ear equations (^^{x^yq) = 0, yg+i) = a].^ {x)y'^_^i. + dial.^{x)y^ = and has symbol 
9q+r = Rq+r ^ Sq-^-r-T* ® E G Jq-f-r{E) if One looks at the top order terms. If /g+i S Rq+i is 
over fq & Rq, differentiating the identity a]!^ {x) f'^{x) = with respect to x^ and substracting the 
identity al^'{x)f^^,^ {x)+dial^'{x)f^{x) = 0, we obtain the identity al>'{x){dj'^(x)-f^+,^ {x)) = 
and thus the restriction D : Rg+i T* ^ Rg ([23], [27], [33]). 

DEFINITION 4.1: Rg is said to be formally integrable when the restriction 7r|+^ : Rg+i — >■ Rg 
is an cpimorphism Vr > or, equivalently, when all the eqiiations of order q + r arc obtained 
by r prolongations only Vr > 0. In that case, Rq+i C Ji{Rq) is a canonical equivalent formally 
integrable first order system on Rg with no zero order equations, called the Spencer form. 

DEFINITION 4.2: Rg is said to be involutive when it is formally integrable and all the sequences 

... A A*T* (gi Qg+r ••• are exact VO < s < n,Vr > 0. Equivalently, using a linear change of lo- 
cal coordinates if necessary, we may successively solve the maximum number (3q , (3q^^ , f3q of 
equations with respect to the principal jet coordinates of strict order q and class n, n — 1, 1 in 
order to introduce the characters = w-f^?rijr(^j3iyT ~ Pg for i = 1, n with — a. Then Rg is 
involutive if -R^+i is obtained by only prolonging the /3* equations of class i with respect to di,...,di 
for 1 = 1, ...,n. In that case dim(gq^i) = + ... + a"^ and one can exhibit the Hilbert polynomial 
dim{Rg+r) in f with leading term {a/n\)r^ when a 7^ 0. Such a prolongation procedure allows 
to compute in a unique way the principal {pri) jets from the parametric {par) other ones. This 
definition may also be applied to nonlinear systems as well. 

We obtain the following theorem generalizing for PD control systems the well known first order 
Kalman form of CD control systems where the derivatives of the input do not appear ([27], VI, 1.14, 
p 802): 

THEOREM 4.3: When Rg is involutive, its Spencer form is involutive and can be modified to a 
reduced Spencer form in such a way that /3 = dim{Rg) — a equations can be solved with respect to 
the jet coordinates 2^, while z^'^^, z^'^" do not appear. In this case z^~^^, z^~^°' do not 
appear in the other equations. 

When Rg is involutive, the linear differential operator "D : E ^ Jq{E) A Jg{E)/Rg = Fq of 
order q with space of solutions © C -B is said to be involutive and one has the canonical linear 
Janet sequence ([4], p 144): 

^ e ^ T A Fo A Fi A ... ^ ^ 

where each other operator is first order involutive and generates the compatibility conditions (CC) 
of the preceding one. As the Janet sequence can be cut at any place, the numbering of the Janet 
bundles has nothing to do with that of the Poincare sequence, contrary to what many physicists 
believe. 

DEFINITION 4.4: The Janet sequence is said to be locally exact at Fr if any local section of F^ 

killed by T>r+i is the image by 'Dr of a local section of Fr-i- It is called locally exact if it is locally 
exact at each Fj. iov Q < r < n. The Poincare sequence is locally exact but counterexemples may 
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exist ([23], p 202). 



Equivalcntly, we have the involutive first Spencer operator Di : Co = Rq ^ Ji{Rq) — > 
Ji{Rq) / Rq+i ^T* ® Rq/5{gq+i) = Ci of Order one induced by D : Rq+i -^T* ® Rg. Introducing 
the Spencer bundles Cr = A^T* ® Rq/d{A'^~^T* gq+i), the first order involutive (r + \)-Spencer 
operator Dr+i : Cr ^ Cr+i is induced by D : N'T* (g) A'^+^T* ® Rq : a Cg+i -> 

da®S,q + {—lya A -D^g+i and we obtain the canonical linear Spencer sequence ([4], p 150): 

^ e ^ Co A Ci C2 ... c„ ^ 

as the Janet sequence for the first order involutive system Rq+\ C J\{Rq). 

The Janet sequence and the Spencer sequence are connected by the following crucial commu- 
tative diagram (1) where the Spencer sequence is induced by the locally exact central horizon- 
tal sequence which is at the same time the Janet sequence for jq and the Spencer sequence for 
Jq+i{E) C Ji{Jq{E)) ([25], p 152): 

SPENCER SEQUENCE 



^ e ^ Co ^ Ci ^ C2 ^...^ c„ -^0 

4- 4- 4' \r 

^ E ^ Co{E) Ci{E) C2{E) Cn{E) ^0 

II i^l i^n 

Q ^ E ^ Fo ^ F, ^ F2 A...^ F„ ^0 
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In this diagram, only depending on ttie left commutative square V = $ o j^, the epimorliisms 
$r : Cr{E) — !> for < r < n arc successively induced by the canonical projection $ = $o : 

Co{E) = Jq{E) ^ Jq{E)/Rq = Fq. 

EXAMPLE 4.5 : {Screw structure): The system i?i C Ji(T) defined by ^2 = 0,^| - S,l = is 
involutive with par{R2) = {^^, Ci ) : Ciu Ciil- The Spencer operator is not involutive as it is 
not even formally integrable because d2£,f — = 0,9i^^ — = ^ ^^iCii ~ '^aCii — 0- We 
obtain dim{Fo) = 2,rfim(Co(T)) = 6 =^ (iim(Co) = dim{Ri) = 4,(iim(Fi) = dim{Ci{T)) = 
dim{Ci) = 6, dim{C2{T)) = dim{C2) = 2 and it is not evident at all that the first order involutive 
operator Di : Cq — s- Ci is defined by the 6 PD equations: 

= 0, d2e - Ci = 0, 92^1^ = 0, 92^1' - d,ei = o, d^e - ei = o, d^e - = o 

The case of a complex structure is similar and left to the reader. 

5 DIFFERENTIAL MODULES AND INVERSE SYSTEMS 

An important but difficult problem in engineering physics is to study how the formal properties 
of a system of order q with n independent variables and m unknowns depend on the parameters 
involved in that system. This is particularly clear in classical control theory where the systems 
are classified into two categories, namely the "controllable" ones and the "uncontrollable" ones 
([14], [27]). In order to understand the problem studied by Macaulay in [M], that is roughly to de- 
termine the minimum number of solutions of a system that must be known in order to determine 
all the others by using derivatives and linear combinations with constant coefficients in a field k, 
let us start with the following motivating example: 
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EXAMPLE 5.1: When n = 1, m = 1, — 3, using ci snb-indcx x for the derivatives with d^y — yx 
and so on, the general solution of yxxx — J/x = is j/ = ae^" + be^^ + cl with a, b, c constants 
and the derivative of is e^, the derivative of e^^ is — and the derivative of 1 is 0. Hence 
we could believe that we need a basis {l,e^ ,e~'^} with three generators for obtaining all the so- 
lutions through derivatives. Also, when n = l,m = 2,k = R and a is a constant real parameter, 
the OD system yl^ — ay^ = 0,yl = needs two generators {(x, 0), (0, 1)} when a = with the 
only dx killing both y^ and 2/2 but only one generator when a ^ 0, namely {{ch{x),l)} when 
a = 1. Indeed, setting y = — y"^ brings y^ = yxxiy"^ = yxx — y and an equivalent system 
defined by the single OD equation yxxx ~ yx = for the only y. Introducing the corresponding 
poynomial ideal (x^ — x) = ix) (x ~ 1) l~l (x + 1), we check that dx kills yxx — y, dx — 1 kills 
Vxx + Vx and dx + l kills yxx~yx, a result leading, as we shall see, to the only generator {ch{x) — 1}. 

More precisely, if ii' is a differential field containing Q with n commuting derivations di, that 
is to say di{a + b) — diO + dib and di{ab) = {dia)b + adib, Va, 6 G for i = 1, n, we denote by 
k a subfield of constants. Let us introduce m differential indeterminates for A; = 1, m and n 
commuting formal derivatives di with diy^ ~ vfi+ii- We introduce the non-commutative ring of 
differential operators D = Kldi, = K[d] with diU — adi + dia,\/a € K in the operator sense 

and the differential module Dy = Dy^ + ... + Dy"\ If {^"^ = a^^y^} is a finite number of elements 
in Dy indexed by r, we may introduce the differential module of equations I = C Dy and the 
finitely generated residual differential module M = Dy/I. 

In the algebraic framework considered, only two possible formal constructions can be obtained 
from M when D = K[d], namely homoiM, D) and M* = homxiM, K) ([3], [27], [32]). 

THEOREM 5.2: homoiM, D) is a right differential module that can be converted to a left dif- 
ferential module by introducing the right differential module structure of a"T* . As a differential 
geometric counterpart, we get the formal adjoint of V, namely ad{V) : A"T* ® F* — ;> a"T* fg) E* 
usually constructed through an integration by parts and where E* is obtained from E by inverting 
the local transition matrices, the simplest example being the way T* is obtained from T. 

REMARK 5.3: Such a result explains why dual objects in physics and engineering are no longer 
tensors but tensor densities, with no reference to any variational calculus. For example the EM 
potential is a section of T* and the EM field is a section of A^T* while the EM induction is a 
section of A^T* (g) A^T :^ A^T* and the EM current is a section of A*T* ® T =i A^T*. 

The filtration Dq = K C Di = K ®T C ... CD, C ... C D of D by the order of operators 
induces a filtration/inductive limit C Mq Q Mi C ... C Mg C ... C M and provides by duality 
over K the projective limit M* = R ^ ... Rg — >■ ... — >■ Ri Rq — >■ of formally integrable 
systems. As D is generated by K and T = Di/Dq, we can define for any / e M*: 

(a/)(m) = a/(m) = /(am), (e/)(m) = ^f{m) - f{^m),ya &K,y^ = a'di &T,^meM 

and check dia = adi + dia, £,r] — ri^ = [£,, 77] in the operator sense by introducing the standard bracket 
of vector fields on T. Finally we get (dif)'^ = {dif){y^) = dif^ — Z^+i^ in a coherent way. 

THEOREM 5.4: R = M* has a structure of differential module induced by the Spencer operator. 

REMARK 5.5: When m = 1 and D = k[d] is a commutative ring isomorphic to the polynomial 
ring ^ = fc[x] for the indeterminates Xii •••iXn, this result exactly describes the inverse system of 
Macaulay with -di = 5i ([M], §59,60). 

DEFINITION 5.6: A simple module is a module having no other proper submodule than 0. A 

semi-simple module is a direct sum of simple modules. When A is a commutative integral domain 
and M a finitely generated module over A, the socle of M is the largest semi-simple submodule 
of M, that is soc{M) = (Bsocm{M) where soCm{M) is the direct sum of all the isotypical simple 
submodules of M isomorphic to A/m for m S max{A) the set of maximal proper ideals of A. The 
radical of a module is the intersection of all its maximum proper submodules. The quotient of a 
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module by its radical is called the top and is a semi-simple module ([3]). 

The " secret " of Macaulay is expressed by the next theorem: 

THEOREM 5.7: Instead of using the socle of M over A, one may use duaUty over k in order to 
deal with the short exact sequence — > rad{R) R ^ top{R) — > where top{R) is the dual of 
soc{M). 

However, Nakayama's lemma ([3], [19], [32]) cannot be used in general unless R is finitely gener- 
ated over k and thus over D. The main idea of Macaulay has been to overcome this difficulty by 
dealing only with unmixed ideals when m = l. As a generalization, one can state ([27]): 

DEFINITION 5.8: One has the purity filtration = tn{M) C ... C to{M) = t{M) C M where 
any involutive system of order p defining Dm is such that a^~'^ = 0, = when m S tr{M) and 
M is said to be r-pure if tr{M) = 0, tr-i{M) = M. With t{M) = {m S M | 30 7^ a £ ^, om = 0} 
we say that M is a 0-pure or torsion-free module if t{M) = and a torsion module if t{M) = M. 

EXAMPLE 5.9: With n = 2,q = 2, let us consider the involutive system y(o,2) = 2/22 = 0, 2/(1,1) = 
J/12 = 0. Then z' = yi satisfies z'2 = while z" = y2 satisfies z'2 = 0, z" = and wc have the 
filtration = t2{M) C ti{M) C ta{M) = t(M) = M with z" G ti{M),z' e to{M) but z' ^ ti{M). 
This classification of observablcs has never been applied to engineering systems like the ones to be 
found in magnetohydrodynamics (MHD) because the mathematics involved are not known. 

REMARK 5.10: A standard result in commutative algebra allows to embed any torsion-free 
module into a free module ([32]). Such a property provides the possibility to parametrize the 
solution space of the corresponding system of OD/PD equations by a finite number of potential 
like arbitrary functions. For this, in order to test the possibility to parametrize a given operator 
T>i, one may construct the adjoint operator ad{'Di) and look for generating CC in the form of an 
operator adiV). As adCD) o adiVi) = adiVi o 2?) = Di o D = 0, it only remains to check 
that the CC of P are generated by Vi. When n = 1 this result amounts to Kalman test and the 
fact that a classical OD control system is controllable if and only if it is parametrizable, a result 
showing that controllability is an intrinsic structural property of a control system, not depending 
on the choice of inputs and outputs contrary to a well established engineering tradition ([14], [27]). 
When n = 2, the formal adjoint of the only CC for the deformation tensor has been used in the 
Introduction in order to parametrize the stress equation by means of the Airy function. This result 
is also valid for the non-commutative ring D = K[d\. 

EXAMPLE 5.11: With K = Q{x^ , x'^ , x^) , infinitesimal contact transformations are defined by 
the system 82^^ - x^d2^^ + x^di^^ - {x^)^di^^ - = 0, 83^^ - x^d^^"^ = 0. Multiplying by test 
functions (A^, A^) and integrating by parts, we obtain the adjoint operator (up to sign): 

d2\^ + x^diX^ + d^X'^ = IJ,\ -x^d2X^ - {x^fdiX^ - x^dsX^ - A^ = /i^ A^ = 

It follows that Ai = fi^, A2 = -/y2 - x^fi^ => d2fi^ + x^difi^ - dsfM^ - x^dsfi^ - 2fj^ = 0. Multi- 
plying again by a test function 0, we discover the parametrization = x^dztj) — 0, = c^af/*, ^'^ = 
—d2(t> — x^dicj) which is not evident at first sight. 

When M is r-pure. Theorem 4.3 provides the exact sequence — M — >■ fc(xi, ...,Xn-r) ® M, 
also discovered by Macaulay ([M], §77, 82), and one obtains the following key result for studying 
the identifiability of OD/PD control systems (see localization in ([19], [27] ,32 [29], [30], [32]). 

THEOREM 5.12: When M is n-pure, one may use the Chinese remainder theorem ([19], p 41) 
in order to prove that the minimum number of generators of R is equal to the maximum num- 
ber of isotypical components that can be found among the various components of soc{M) or top{R). 
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6 JANET VERSUS SPENCER : THE NONLINEAR SEQUENCES 



Nonlinear operators do not in general admit CC as can be seen by considering the involutive 
example 2/22 — 5(2/11)"^ = u, yi2 — \{yii)^ = v with m = 1, n = 2, g = 2, contrary to what happens 
in the study of Lie pseudogroups. However, the kernel of a linear operator V : E ^ F is always 
taken with respet to the zero section of F, while it must be takc;n with respect to a prescribed 
section by a double arrow for a nonlinear operator. Keeping in mind the linear Janet sequence and 
the examples of Vessiot structure equations already presented, one obtains: 

THEOREM 6.1: There exists a nonlinear Janet sequence associated with the Lie form of an 
involutive system of finite Lie equations: 

*w Ojq I O ji 

-)■ r aui(X) T ^ J^i 

LOO a 

where the kernel of the first operator / — !■ $0; ° jqif) = ^Lj{jq{.f)) = jgi.f)"^ {^) is taken with re- 
spect to the section uj of while the kernel of the second operator w — )■ I{ji{ui)) = A{ui)dx(jJ + B{ui) 
is taken with respect to the zero section of the vector bundle J^i over J^. 

COROLLARY 6.2: By linearization at the identity, one obtains the involutive Lie operator 
D = D„ : T Fo : ^ -)■ C{^)co with kernel 6 = G T\C{^)w = 0} C T satisfying [6, 6] C 6 and 
the corresponding linear Janet sequence where Fq = uj~^{V{J^)) and Fi = oj~^{J^i). 

Now we notice that T is a natural vector bundle of order 1 and Jq{T) is thus a natural vector 
bundle of order q+ 1. Looking at the way a vector field and its derivatives are transformed under 
any / G aut{X) while replacing jq{f) by fg, we obtain: 

and so on, a result leading to: 

LEMMA 6.3: JqiT) is associated with Ilg+i = Ilq+i{X, X) that is we can obtain a new section 
r]q = fq+i{^q) from any section ^q G Jq{T) and any section G Hq+i by the formula: 

d^ri'' = r^^fl + ... = + ... + 4\i,r,V0 < ImI < 9 
where the left member belongs to ^(ng). Similarly Rq c Jg(T') is associated with 7^g+i C 11,+ 1. 

In order to construct another nonlinear sequence, we need a few basic definitions on Lie 
groupoids and Lie algebroids that will become substitutes for Lie groups and Lie algebras. As in the 
beginning of section 3, the first idea is to use the chain rule for derivatives jq{g o /) = jq{g) o jq{f) 
whenever /, 5 G aut{X) can be composed and to replace both jq{f) and jq{g) respectively by fq 
and gq in order to obtain the new section gq o fq. This kind of "composition" law can be written in 
a pointwise symbolic way by introducing another copy Z of X with local coordinates (z) as follows: 

jq : nq{Y, Z)xynqiX, Y) ^ n,(X, Z) {{y, z, g, ...), {X, y, g, ...) ^ (a;, ^, |^ ...) 
We may also define jq{f )~^ = jq{f~^) and obtain similarly an "inversion" law. 

DEFINITION 6.4: A fibered submanifold TZq C Ilg is called a system, of finite Lie equations 
or a Lie groupoid of order q if we have an induced source projection Oq : TZq ^ X , target pro- 
jection j3q : TZq X, composition 7, : TZqXxTiq Ttq, inversion Lq : TZq TZq and identity 
idq : X ^ TZq. In the sequel we shall only consider transitive Lie groupoids such that the map 
{aq,f3q) : Kq ^ X X X is an epimorphism and we shall denote by TZq = id~^{TZq) the isotropy 
Lie group bundle of TZq. Also, one can prove that the new system priJZq) = T^q+r obtained by 
differentiating r times all the defining equations of TZq is a Lie groupoid of order q-\- r. Finally, 
one can write down the Lie form and obtain TZq = {fq G 11,1/^^(0;) — w}. 
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Now, using the algebraic bracket {jq+iiO^ Jq+iiv)} = jq{[^jV])i^^''l S T, we may obtain by 
bilinearity a differential bracket on Jq{T) extending the bracket on T: 

[(.Q-Vq] = {^q+l,Vq+l} + KO^Vq+l " «(»7)^^9+l , V^g , ??q G Jq{T) 

which does not depend on the respective Ufts ^g+i and r]q+i of and rjq in Jq+i(T). This bracket 
on sections satisfies the Jacobi identity and we set: 

DEFINITION 6.5: We say that a vector subbundle Rq C Jq{T) is a system of infinitesimal Lie 
equations or a Lie algebroid if [Rq,Rq] C Rq, that is to say [^q,r]q] G Rq,y^q,r]q G Rq. The kernel 
iig of the projection ttq : Rq ^ T is the isotropy Lie algebra bundle of 7?.° and R^] C ii^ does 
not contain derivatives. 

PROPOSITION 6.6 : There is a nonlinear diff'erential sequence: 

— > aut{X) ng+i{X,X) A T* (g) Jg{T) ^ A^T* ® J,_i(T) 

withl)/g+i EE fq^ioji{fg)-idq+i =Xq^ £>'xq{^,v) = DxqH, v) - {XqiO , Xqiv)} = 0- Moreover, 
setting xo = A — id £T* ^T, this sequence is locally exact if det{A) ^ 0. 

Proof: There is a canonical inclusion llg+i C Ji(i^q) defined by y^^^ = y^+i- and the composition 
fq+i o jiifq) is a well defined section of JiiTlq) over the section f~^ ° fq = idq of Ilg like idq+i. 
The difl:erence Xq = fq+i ° Ji(/g) ~ idq+i is thus a section of T* ® X''(n,,) over id, and we have 
already noticed that id~^{V{Hq)) = Jq{T). For g = 1 we get with gx = /f^: 

X% = 9!tdif - = A\ - 5l xli = gHdifj - A^lj) 

We also obtain from Lemma 6.3 the useful formula frXli,i + ••• + f^+i^X^i — ^if^ ~ f^+u allowing 
to determine Xq inductively. 

We refer to ([26], p 215) for the inductive proof of the local exactness, providing the only formulas 
that will be used later on and can be checked directly by the reader: 

^iX,j ~ ^jX^i ^ XiJ ^" Xj,i ^ iX.iXr.j ^ X.jXr.i) ~ 
~ ^jXl.i XliJ + Xlj.i {X,iXlr,j ~^ Xl.iXrJ ^ Xl.jXr.i ^ XjXlr.i) ^ 

There is no need for double-arrows in this framework as the kernels are taken with respect to the 
zero section of the vector bundles involved. We finally notice that the main difference with the 

gauge sequence is that all the indices range from 1 to n and that the condition det{A) amounts 
to A = det{dif'^) ^ because det{f^) 7^ by assumption. □ 

COROLLARY 6.7: There is a restricted nonlinear diff'erential sequence: 

^ r Tlq+l A T* ® i?, A^T* Jq-l{T) 

DEFINITION 6.8: A splitting of the short exact sequence ii° -> i?, T is a map x'q ■ 
T ^ Rq such that ttq 0%^ = idr or equivalently a section of T* ®Rq over idr gT*^T and is called 

a Rq- connection. Its curvature k'^ e A^T* ® i?" is defined by K'q{(,,r]) = [x'qiO-x'qiv)] ~ Xg([^)^])- 
We notice that x!q = —Xq is a connection with D'xlq = if and only if ^ = 0. In particular 
(5f , is the only existing symmetric connection for the Killing system. 

REMARK 6.9: Rewriting the previous formulas with A instead of xo we get: 

d,A) - d^A\ - A\x\. + A^r^K = 

^ixt- - ^jXt - XmX^j + xLXr.i - ^Ixfr,,- + ^Jxt.i = 

When q= \,gi = and though surprising it may look like, we find back exactly all the formulas 
presented by E. and F. Cosserat in ([C], p 123 and [16]). Even more strikingly, in the case of a 
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Riemann structure, the last two terms disappear but the quadratic terms are left while, in the case 
of screw and complex structures, the quadratic terms disappear but the last two terms are left. 

COROLLARY 6.10: When det{A) ^ there is a nonlinear stabilized sequence at order q: 

aut{X) ^ n, Ci(T) C2{T) 
called nonlinear Spencer sequence where Di and D2 are involutive and its restriction: 

r ^ 7^g A Ci A C2 

is such that Di and D2 are involutive whenever TZg is involutive. 

Proof: With = g we have ^ = — gf fj^^i^+terms{order < q). Setting ^ = A^t^ .,., we ob- 
tain r^,. = -gf f]j^^-i^^+terms{order < q) and D : n,+i -y T*®Jq{T) restrics to 5i : n, C'i(r). 
Finally, setting A~^ = B = id — tq, we obtain successively: 

^^xiJ - ^Jxi^ + termsixg) ~ (A[x^+i,,j - A^jXl+u,r) = 
BlBl{dixl,j - ^,xi^) + termsixq) - (t^+i,,. - r^+i^^r) = 

We obtain therefore Dr^+i + terms{Tg) = and D' : T* ® JqiT) A^T* (gi Jq-i{T) restricts to 
D2 : CiiT) ^ C2iT). 

In the case of Lie groups of transformations, the symbol of the involutive system Rq must be = 
providing an isomorphism TZq+i ~ 7^^ Rq+i — Rq and we have therefore Cr = A'^T* Rq for 
r = l,...,n in the linear Spencer sequence. □ 

REMARK 6.11: The passage from Xg to r, is exactly the one done by E. and F. Cosserat in ([C], 
p 190). However, even if is a good idea to pass from the source to the target, the way they realize 
it is based on a subtle misunderstanding that we shall correct later on in Proposition 6.16. 

If fq+i,9q+i e n,+i and = o wc get: 

£>f'q+l = fq'+l ° 9gll ° Jl (ffg) ° jl ifq) " idq+1 = fg+l O Dgg+l O ji (/,) + £»/,+ ! 

DEFINITION 6.12: For any section fg+i € TZg+i, the transformation: 

Xq >Xq = fg+l °Xq° hifq) + Dfq+1 

is called a gauge transformation and exchanges the solutions of the field equations D'xg = 0. 
Introducing the formal Lie derivative on Jg(T) by the formulas: 

L{^q+l)'nq = {iq+l,Vq+l} + iiO^nq+l = [^q, Vq] + i{'n)D^.q+l 
iL{h{^q+l))Xqm = ife+l)(X,(C)) -X,(K,C]) 

and passing to the limit with /g+i = idq+i + t^q+i + ... for t — >^ over the source, we get: 
LEMMA 6.13: An infinitesimal gauge transformation has the form: 

SXq = + L{jl{^q+l))Xq 

Passing again to the limit but now over the target with Xq = -D/g+i and t/g+i = idq+i+tr]q+i+..., 
we obtain the variation: 

SXq = .fq'+l O Dr]q+l O ji{fq) 

PROPOSITION 6.14: The same variation is obtained whenever rjg+i — fg+2{£,g+i + Xg+i(C)) 
with Xq+i = Dfq+2, a transformation which only depends on ji{fq+i) and is invertible if and only 
if det{A) ^ 0. 
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Proof: Choosing /g+i , /i^+i £ Tlq+i such that g^+i o f^^i = f^^i o /i^+i and passing to the 
limits Qq+i = idq+i + triq+i + ... and hq+i — idq+i + + ■•■ when i 0, we obtain the local 
formula: 

^M^' = ^'/; + - = e{^^fl - + + - + fUi 

and thus r/g+i = /5+2(f9+i) with ^g+i — ^g+i + Xg+i(C)- This transformation is invertible if and 
only if ^ — )• ^ = ^ + ;)(;o(^) = A(^) is an isomorphism of T. □ 

EXAMPLE 6.15: For g = 1, we otain: 

<5xfi = {die - e) + iedrx'^i + x':Ac - x:ier) 

Sxii = - ^j) + i^drxl^ + xlAC + X.'/I -Xli^r - x:i^i) 

i^i^j ^ij) ~^ iXrj.i^ ~^ Xr.i^j Xj,i^r X.i^jr^ 

For the Killing system Ri C Ji{T) with g2 = 0, these variations are exactly the ones that can 
be found in ([C], (50)+(49), p 124 with a printing mistake corrected on p 128) when replacing a 
3x3 skewsymmetric matrix by the corresponding vector. The last unavoidable Proposition is thus 
essential in order to bring back the nonlinear framework of finite elasticity to the linear framewok 
of infinitesimal elasticity that only depends on the linear Spencer operator. 
For the conformal Killing system .Ri C Ji (T) (see next section) we obtain: 

This is exactly the variation obtained by Weyl ([W], (76), p 289) who was assuming implicitly 
A = when setting =0 ^ = — a^^* by introducing a connection. Accordingly, is the vari- 
ation of the EM potential itself, that is the 8Ai of engineers used in order to exhibit the Maxwell 
equations from a variational principle ([W], § 26) but the introduction of the Spencer operator is 
new in this framework. 

Finally, chasing in diagram (1) , the Spencer sequence is locally exact at C\ if and only if 
the Janet sequence is locally exact at Fq because the central sequence is locally exact. The sit- 
uation is much more complicate in the nonlinear framewok. Let w be a section of T satisfying 
the same CC as w, namely /()i(w)) = 0. It follows that we can find a section fq+i £ Hg+i such 
that f-\u:) Jiifq'){ji{io)) = JiifqH^)) = 3i{0j) and /-+\(ji(a;)) = ji(w). We obtain 

therefore ji_(/-i)(ji(w)) = /-+\(ji(w)) ^ Uq+i ° Ji{fq^))-\ji{^)) - 3i{^) = L{aq)uj = and 
thus cTg = Df~^-^ e T* i?g over the target, even if /g+i may not be a section of T^^+i. As a, is 
killed by D', wc have related cocycles at in the Janet sequence over the source with cocycles at 
T* (g) Rq or Ci over the target. 

Now, if /,+i,/g_^i G Ilg+i are such that /g+\(ji(w)) = fq+\{ji{u))) = ji(w), it follows that 
i.fq+i°.fq+i)i.iii^^)) = .?i('^) => 3.9g+i G Uq+i such that =5(,+io/g+i and the new ct^ = Df'g~\ 
differs from the initial aq = Df^^^ by a gauge transformation. 

Conversely, let fq+i,fq_^_i € ^q+i be such that aq = Df~_^^ = Df'^ ~ '^r follows that 
D{fq+x o f'^^^) = and one can find g G aut{X) such that /g+i = fq+i o jq+i{g) providing 

U>' = f'q-^u) = {fqOjq{g))-\uj)=jq{gr\fq-\uj))=jq{g)-\uj). 

PROPOSITION 6.16: Natural transformations of over the source in the nonlinear Janet 
sequence correspond to gauge transformations of T* (g) Rq or Ci over the target in the nonlin- 
ear Spencer sequence. Similarly, the Lie derivative = C{£,)io G Fq in the linear Janet sequence 
corresponds to the Spencer operator -D^g+i G T*^Rq or Di^q G Ci in the linear Spencer sequence. 

With a slight abuse of language Sf = r]of <^ Sfof~^ = rj f^^°5f = ^ when r) = T[f){^) and 
weget j,(/)-i(a;) = uj ^ jq{f +5f)-^{uj) = Q+Suj fh&tis 3q{f-^o{f +5f))-\Lo) = w-l-to ^ = 
£(0^andjg((/ + (5/)o/-io/)-i(a;)=i,(/)-i(j,((/ + 5/)o/-i)-i(a;))^<5w=j,(/)-i(£(r?)a;). 
Passing to the infinitesimal point of view, we obtain the following generalization of Remark 3.12 
which is important for applications ([2], AJSE-mathematics): 
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COROLLARY 6.17: 5u) = C{^)uj = j,(/)-i(£(ry)w). 



EXAMPLE 6.18: In Example 3.4 with n = l,q = 1, we have w(/(a;))/,;(x) = u}{x),uj{f{x))f^^{x) + 
dyLo{J{x))fl{x) = dxUj{x) and obtain therefore oJCJy^y + a^ydyU! = -uj{l/ fx){dxfx - fxx){^/dxf) + 
{{.fx/dxf) — l)dyUJ = whenever y = f{x). The case of an afHne stucture needs more work. 

7 COSSERAT VERSUS WEYL: NEW PERSPECTIVES FOR PHYSICS 

As an application of the previous mehods, let us now consider the conformal Killing system: 

Rl C Jl{T) U)rj(.l + tOirCj + Cdrl^ij = A.{x)'^i3 

with symbols: 

obtained by eliminating the arbitrary function A{x), where ui is the Euclidean metric when n = 1 
(line), n = 2 (plane) or n = 3 (space) and the Minskowskian metric when n = 4 (space-time). 
The brothers Cosserat were only dealing with the Killing subsystem: 

Rl C -Rl COrj^l + OOir^j + Cdr^ij = 

that is with Ci I = ^'^ij = 0} = {translations, rotations} when A{x) = 0, while, in a 
somehow complementary way, Weyl was mainly dealing with {^^'Ci} = {dilatation, elations}. 
Accordingly, one has ([7]): 

THEOREM 7.1: The Cosserat equations ([C], p 137 for n = 3, p 167 for n = 4): 

are exactly described by the formal adjoint of the first Spencer operator Di : Ri ^ T* ® Ri. 
Introducing (jf^^^ = —<j)^'^'^ and tp'^^''^^ = — = —ip^^^^i^ they can be parametrized by the formal 
adjoint of the second Spencer operator D2 : T* ® Ri ^ A^T* Ri: 



EXAMPLE 7.2: When n = 2, lowering the indices by means of the constant metric u), we just 
need to look for the factors of ^1,^2 and ^1,2 in the integration by parts of the sum: 

(T''\di^l - a,l) + CT''\d2^1 - a,2) + CT^'HdiC2 ' + T^'\d2^2 - 6,2) + fJ-'^'"' {drCl,2 - Cl,2r) 

Finally, setting 0^^^^ = — — <j)^ , we obtain the nontrivial parametrization 

al.l = S2<^l,al'2 = -5i0l,(72'l = -92.^2^(72-2 = ai.^2^^12,1 ^ Q^^Z ^ ^i^^-^ = -d^j}^ - (/.2 in a 

coherent way with the Airy parametrization obtained when (f>^ = d2(f>, (fy^ = di(f>, cffi = —(f>. 

REMARK 7.3 : First of all, it is clear that [C] (p 13,14 for n = 1, p 75,76 for n = 2) still deals 
with m = 3 for the "ambient space", that is with the construction of the nonlinear gauge sequence, 
in particular for the dynamical study of a line with arc length s and time t considered as a surface, 
hence with no way to pass from the source to the target, only possible, as we have seen, when 
TO = n = 3 by using the nonlinear Spencer sequence. For n = 4, the group of rigid motions of space 
is extended to space-time by using only a translation on time and we can rewrite the formulas in 
([C], p 167) as follows: 

d dx d dy d dz d d dpxx 1 dA d A dx /'^\ 

di " 'dtd^ did^ dtd^ ^dt^ ^ ■■■ A " d^^^'"' ^ A ^ Wa' 

dQxx IdP Cdy Bdz d, Pdx, d,P^ , Cdy, , Bdz, 

-^+-+Py^-P^y+A^+Adi-Adi = a^^«^^+A^^+-+at(A)+(^-+A rf^)-(^-+A^) 
It is essential to notice that the Cosserat equations for n = 3 are still introduced today in a phe- 
nomenological way ([35] is a good example), contrary to the "deductive" way used in ([C], p 1-6) 



20 



and that "intuition" will never allow to provide the relativistic Cosserat equations for n = 4 which 

arc presented for the first time. 

THEOREM 7.4: The Weyl equations ([W], §35) are exactly described by the formal adjoint of 
the first Spencer operator Di : R2 ^ T* ^ R2 when n = 4 and can be parametrized by the for- 
mal adjoint of the second Spencer operator D2 : T* (g) R2 ^ T* (E) R2- In particular, among the 
components of the first Spencer operator, one has di£,^j — ^[^^ = Sj^^^- and thus the components 
di£,rj — dj£,^i = Fij of the EM field with EM potential = Ai coming from the second order jets 
(elations). It follows that Di projects onto d : T* ^ A^T* and thus D2 projects onto the first 
set of Maxwell equations described by d : A^T* A^T*. Indeed, the Spencer sequence projects 
onto the Poincarc sequence with a shift by +1 in the degree of the exterior forms involved because 
both sequences are made with first order involutivc operators and the comment after diagram (1) 
can thus be used. By duality, the second set of Maxwell equations thus appears among the Weyl 
equations which project onto the Cosserat equations because of the inclusion i?i ~ i?2 C i?2- 

REMARK 7.5: When n — A, the Poincare group (10 parameters) is a subgroup of the conformal 
group (15 parameters) which is not a maximal subgroup because it is a subgroup of the Weyl 
group (11 parameters) obtained by adding the only dilatation with infinitesimal generators x^di. 
However, the optical group is another siibgroup with 10 parameters which is maximal and the same 
procedure may be applied to all these subgroups in order to study coupling phenomena. It is also 
important to notice that the first and second sets of Maxwell equations are invariant by any difjeo- 
morphism and the conformM group is only the group of invariance of the Minkowski constitutive 
laws in vacuum ([20]) ([27], p 492). 

REMARK 7.6: Though striking it may look like, there is no conceptual difference between the 
Cosserat and Maxwell equations on space-time. As a byproduct, separating space from time, there 
is no conceptual difference between the Lame constants (mass per unit volume) of elasticity and 
the magnetic (dielectric) constants of EM appearing in the respective wave speeds. For example, 
the speed of longitudinal free vibrations of a thin elastic bar with Young modulus E and mass 

per unit volume p is v = while the speed of light in a medium with magnetic constant fj, and 
dielectric constant e is w = \J~^^- In the first case, we have the 1-dimensional dynamical equations: 

In the second case, studying the propagation in vacuum for simplicity, one uses to set H = 
{l/fiQ)B,D = cqE with eonoc^ = 1 in the induction equations and to substitute the space-time 
parametrization dA = F oi the field equations dF = in the variational condition S J{^eoE'^ — 
^{l/lJ,o)B'^)dxdt = 0. However, the second order PD equations thus obtained become wave equa- 
tions only if one assumes the Lorentz condition div{A) = uj^^diAj = ([20]). This is not correct 
because the Lagrangian of the corresponding variational problem with constraint must contain the 
additional term \div{A) where A is a Lagrange multiplier providing the equations DA = dA as a 
1-form and thus UF = as a 2-form when □ is the Dalembertian ([27], p 885). 

REMARK 7.7: When studying static phenomena, e = {cij) and E = [E^) are now on equal 
footing in the Lagrangian, exactly like in the technique of finite elements. Starting with a homoge- 
neous medium at rest with no stress and electric induction, we may consider a quadratic Lagrangian 
A^^^^CijCki+B'^^ EiEj + C^^^tijEk obtained by moving the indices by means of the Euclidean metric. 
The two first terms describe (pure) linear elasticity and electrostatic while only the last quadratic 
coupling term may be used in order to describe coupling phenomena. For an isotropic medium, 
the 3-tensor C must vanish and such a coupling phenomenon, called piezzoelectricity, can only 
appear in non-isotropic media like crystals, providing the additional stress a*^ = C^^^Ek and/or 
an additional electric induction — C^^^eij. Accordingly, if the medium is fixed, for example 
between the plates of a condenser, an electric field may provide stress inside while, if the medium is 
deformed as in the piezzo-lighters, an electric induction may appear and produce a spark. Finally, 
for an isotropic medium, we can only add a cubic coupling term C^^'^'-CijEkEi responsible for pho- 
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toelasticity as it provides the additional electric induction = [C^^^^ eij)Ek, modifying therefore 
the dielectric constant by a term depending linearly on the deformation and thus modifying the 
index of refraction n because e/ioc^ = v? with eoMoc^ = 1 in vacuum leads to e = n^eo. We may 
also identify the dimensionless "speed" w'^/c ^ 1, VA: = 1, 2, 3 (time derivative of position) with a 
first jet (Lorentz rotation) by setting 84^'^ — ^| = and introduce the speed of deformation by the 

formula 2l^y = UJrj{di^l - ^'4) + ^^ir{djS,l - CJ4) = iOrjdi^l + Uirdj^l = di{uJrjdiC + ^irdjC) = 

LOrjdiV^ + LOirdjv"^ = d^Cij < 1, j < 3 in order to obtain streaming birefringence in a similar way. 
These results perfectly agree with most of the field-matter couplings known in engineering sciences 
([28]) but contradict gauge theory ([15], [26]) and general relativity ([W],[21]). 

In order to justify the last remark, let G be a Lie group with identity e and parameters a 
acting on X through the group action X x G ^ X : (x^a) ^ y = f{x, a) and (local) infinitesimal 
generators 9r satisfying [9p,6a] = c^^9t for p,a,T = 1, ...,dim{G). We may prolong the graph 
of this action by differentiating q times the action law in order to eliminate the parameters in 
the following commutative and exact diagram where TZq is a Lie groupoid with local coordinates 
{x,yq), source projection aq : {x,yq) — {x) and targiei projection /?, : {x,yq) — >■ (y) when q is large 
enough: 

0^ X xG — > Tlq -)-0 

II / \ /3g 

X XG ^ X X X 

The link between the various sections of the trivial principal bundle on the left {gauging proce- 
dure) and the various corresponding sections of the Lie groupoid on the right with respect to the 
source projection is expressed by the next commutative and exact diagram: 

Oh- X X G = Uq -^0 

a = cstmaix) jqif)Utfq 
X = X 

THEOREM 7.8: In the above situation, the nonlinear Spencer sequence is isomorphic to the 
nonlinear gauge sequence and we have the following commutative and locally exact diagram: 

X xG ^ T*®g ^ A^T* (g) g 

4-4- 4^ 

O-^r^ Uq ^ T*®Rq ^ A'^T*(E)Rq 

The action is essential in the Spencer sequence but disappears in the gauge sequence. 

Proof: If we consider the action y = f{x,a) and start with a section (x) — >■ {x,a{x)) of X x G, 
we obtain the section {x) — >■ {x,f^{x) = dij,f^{x,a{x))) of TZq. Setting b = = b{a), we get 
y = fix, a) ^x = f{y,b) ^y = /(/(y, 6(a), a) and thus ||§£||+ |m = o with f = e{x)w{b) 
from the first fundamental theorem of Lie. With —co{b)db = —dbb ^ = a ^da, we obtain: 

dif^ - f^+u = diid^^fHx, a{x)) - d^+ufix, a{x)) 

= -d,{^ei{x))u^l{b)^d,a^ 

and thus X^^i{x) = AJ {x)dn9^{x) from the inductive formula allowing to define Xq = ^fq+i- 
As for the commutatitvity of the right square, we have: 

({x,+i(a,),x,+i(a,)})^ = KA]d^{\9,,9,\f = c;,AfAJ5^0^ □ 

Introducing now the Lie algebra Q = Te{G) and the Lie algebroid Rq C Jq{T), namely the 
linearization of TZq at the g'-jet of the identity y = x, we get the commutative and exact diagram: 
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0^ X X g = Rg -^0 

x = cstmm ig(ontc« 

X = X 

where the upper isomorphism is described by X'^{x) C^(a^) = {x)d^9'^{x) for q large enough. 
The unusual Lie algebroid structure on X x ^ is described by the formula: ([A, A'])^ ~ c^^A^A'"^ + 
{\P6p).\''^ — {\''^d„).X^ which is induced by the ordinary bracket on T and thus depends on 

the action. Applying the Spencer operator, we finally obtain diS,^,{x) — (x) = diX'^ {x)dij,9^{x) 
and the linear Spencer sequence is isomorphic to the linear gauge sequence already introduced 
which is no longer depending on the action as it is only the tensor product of the Poincare se- 
quence by g. 

EXAMPLE 7.9: Let us consider the group of affine transformations of the real line y = a^x + a} 
with n = \,dira(G) = 2,q = 2, Tl2 defined by the system y^x = 0, R2 defined by S^^x = and 
the two infinitesimal generators 9i = x-^,92 = -§^- We get f{x) — a^{x)x + a^{x), fx{x) = 
fxx{x) = and thus x x{x) = {1/ fx{x))dxf {x) — 1 = {l/a^{x)){xdxa^{x) + dxO.'^ix)) = 
xAlix) + Al{x),XxAx) = {^/Ux)){dxfx{x) - {l/fx{x))dxf{x)fxx{x)) = {l/a\x))dxa\x) = 
Al{x),Xxx,x{x) = 0. Similarly, we get ^{x) = X^{x)x + X'^{x),^x{x) = X'^{x),^xx{x) = 0. Finally, 
integrating by part the sum (T{dxS, — £,x) + IJ-{dxCx — £,xx) we obtain the dual of the Spencer opera- 
tor as dxO- = /, dxH+cr = m that is to say the Cosserat equations for the affine group of the real line. 

It finally remains to study GR within this framework, as it is only " added" by Weyl in an 
independent way and, for simplicity, we shall restrict to the linearized aspect. First of all, it be- 
comes clear from diagram (1) that the mathematical foundation of GR is based on a confusion 
between the operator Vi [classical curvature alone) in the Janet sequence when V is the Killing 
operator brought to involution and the operator D2 {gauge curvature=curvature+torsion) in the 
corresponding Spencer sequence. It must also be noticed that, according to the same diagram, 
the bigger is the underlying group, the bigger are the Spencer bundles while, on the contrary, the 
smaller are the Janet bundles depending on the invariants of the group action (deformation tensor 
in classical elasticity is a good example). Precisely, as already noticed in Theorem 7, if G C G, the 
Spencer sequence for G is contained into the Spencer seqiience for G while the Janet sequence for 
G projects onto the Janet sequence for G, the best picture for understanding such a phenomenon 
is that of two children sitting on the ends of a beam and playing at see-saw. 

Such a confusion is also combined with another one well described in ([40], p 631) by the Chinese 
saying " To put Chang 's cap on Li 's head'' , namely to relate the Ricci tensor (usually obtained from 
the Riemann tensor by contraction of indices) to the energy- momentum tensor (space-time stress), 
without taking into account the previous confusion relating the gauge curvature to rotations only 
while the (classical and Gosserat) stress has only to do with translations. In addition, it must be 
noticed that the Cosserat and Maxwell equations can be parametrized while the Einstein equations 
cannot be parametrized ([29]). 

In order to escape from this dilemna, let us denote by B'^{gq), Z'^(gq) and H'^{gq) = Z'^(gq)/B'^{gq) 
the coboundary (image of the left S), cocycle (kernel of the right 6) and cohomology bundles of the 

^-sequence T* (g) gq+i A A^T* ® 5^ A A^T* (g) Sq-iT* (g) T. It can be proved that the order of the 
generating CC of a formally integrable operator of order q is equal to s -I- 1 when s is the smallest 
integer such that H'^{gq+r) = 0,Vr > s ([26]). As an example with n = 3, we let the reader prove 
that the second order systems 1/33 = 0,2/23-2/11 = 0,2/22 = and 2/33-2/11 = 0,2/23 = 0,2/22-2/11 = 
have both three second order generating CC ([30]). For the Killing system Ri C Ji{T) with symbol 
gi, we have Fq = Ji{T)/Ri = T*®T/g-i and the short exact sequence ->• gi T*®T -)• 0. 
As = 1 and 52 = 33 = we have 8 = 1 and no GG of order 1. The generating GG of order 
2 only depend on Fi = w~^(J^i) according to section 2 where Fi is now defined by the following 
commutative diagram with exact columns but the first on the left and exact rows: 
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It follows from a chasc([26], p 55)([27], p 192)([32], p 171) that there is a short exact connect- 
ing sequence — > B^{gi) — >• Z^{gi) ^ Fi — >• leading to an isomorphism Fi ~ H^{gi). The 
Riemann tensor is thus a section of Riemann = Fi = H'^{g\) = Z'^{g\) in the Killing case with 
dim{Riemann) = (n2(n + 1)2/4) - (^^(n + l)(n + 2)/6) = (^^(n - 1)^/4) - (^^(n - l)(r7, - 2)/6) = 
n^{n^ — 1)/12 by using either the upper row or the left column and we find back the two algebraic 
properties of the Riemann tensor without using indices. 

However, for the conformal Killing system, we still have q = 1 but the situation is much more 
delicate because (?3 = for n > 3 and H'^{g2) = only for n > 4 (]26], p 435). Hence, setting 
similarly Fq = T* T/51, the Weyl tensor is a section of Weyl = Fi = H'^{gi) 7^ Z'^{gi). The 
inclusion gi C gi and the relations 32 = 0, 53 = finally induce the following crucial commutative 
and exact diagram (2) ([25], p 430): 





Ricci 

1 ; 

Z'^{gi) ^ Riemann -J> 

; i i JANET 

T*®g2 4 Z2(^i) ^ Weyl ^0 

4- •J' 4- 

S2T* ^ T*®T* 4 A^T* -> 

1 i 


SPENCER 



A diagonal chase allows to identify Ricci with S2T* without contracting indices and provides 
the splitting of T* (g T* into 5*27^* {gravitation) and A^T* {electromagnetism) in the lower hori- 
zontal sequence obtained by using the Spencer sequence, solving thus an old conjecture. However, 
T* T* c:i T* (E) g2 has only to do with second order jets {elations) and not a word is left from the 
standard approach to GR. In addition, we obtain the following important theorem explaining for 
the first time classical results in an intrinsic way: 

THEOREM 7.10: There exist canonical splittings of the various i5-maps appearing in the above 
diagram which allow to split the vertical short exact sequence on the right. 

Proof: We recall first that a short exact sequence M' 4 M A M" of modules splits, 
that is M ~ M' ® M" , if and only if there exists a map u : M ^ M' with u o f = idu' or a 
map w : M" ^ M with gov = idM" ([3] p 73) ([32], p 33). Hence, starting with (t,^^.) £ T* ig ^2, 
we may introduce {pf ^j = rj^j - rj'j ^) £ B'^{gi) C Z'^{gi) C A^T* ® gi but now = p\.^j = 
- Trj.i = Pij - Pji 7^ with pij = plj,j because we have p^'j^ + p*",^, + p^j^ = 0. With 
r = w'^T^j J and p = cj'^^pij, we obtain (n — 2)t^^j = (n — l)pij + pji — (n/2(n — l))iOijp and thus 
np = 2(n — 1)t. The lower sequence splits with (pij — >• Tij = t^^j = {l/2)ipij — >• Ty — Tji = ipij and 
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Pij — Pji fij ~ ill ^^(.9i) ^ A^T* It follows from a chase that the kernel of the canonical 
projection Riemann — ^ Weyl is defined by p^j = t^^^ — Ty ^ with (pf j^) G .^^((71) C Z^{gi) and 
(tj^j) G T* (g) ^2- Accordingly (n — 2)rij = npij — (n/2(n — provides the isomorphism 

Ricci ~ S2T* and we get npf ^^- = 5fry - ^j^r;, + wijU^^Tgi - uJuw^^Tsj, that is: 

pin = -^^^r^^^^PiJ - ^jP^^ + ^'J^'^'P- - ^u^'^'Psj) - (^n-l){n-2) ^^'^'^'' ~ ^^"^''^^ 

We check that p[^^^- = pij , obtaining therefore a splitting of the right vertical sequence in the last 
diagram that allows to define the Weyl tensor by difference. These purely algebraic results only 
depend on co independently of any conformal factor. □ 

EXAMPLE 7.11: The free movement of a body in a constant static gravitational field g is de- 
scribed by ^ ^ = 0, ^ — g = 0, ^ --0 = where the "speed" is considered as a first order 
jet (Lorcntz rotation) and the "gravity" as a second order jet (elation). Hence an accelerometer 
merely helps measuring the part of the Spencer operator dealing with second order jets {equivalence 
principle). As a byproduct, the difference 94/4 — /I4 under the constraint 84/'' — identifying 
the "speed" with a first order jet allows to provide a modern version of the Gauss principle of 
least constraint where the extremum is now obtained with respect to the second order jets and not 
with respect to the "acceleration" as usual ([1], p 470). The corresponding infinitesimal variational 
principle 5 J (p(94^^ — ^4) + g^{diSJ^ — ^^i) + 9^'' idiS,rj — Q))dx = provides the Poisson law of gravi- 
tation with p = est and g = (g*) when g*-' = Aw'-' =^ Qi = —diX. The last term of this gravitational 
action in vacuum is thus of the form Xdiv{A), that is exactly the term responsible for the Lorentz 
constraint in Remark 7.6. 

8 CONCLUSION 

In continuum mechanics, the classical approach is based on differential invariants and only 
involves derivatives of finite transformations. Accordingly, the corresponding variational calculus 
can only describe forces as it only involves translations. It has been the idea of E. and F. Cosserat 
to change drastically this point of view by considering a new differential geometric tool, now called 
Spencer sequence, and a corresponding variational calculus involving both translations and rota- 
tions in order to describe torsors, that is both forces and couples. 

About at the same time, H. Weyl tried to describe electromagnetism and gravitation by using, 
in a similar but complementary way, the dilatation and elations of the conformal group of space- 
time. We have shown that the underlying Spencer sequence has additional terms, not known today, 
wich explain in a unique way all the above results and the resulting field-matter couplings. 

In gauge theory, the structure of electromagnetism is coming from the unitary group U{1), 
the unit circle in the complex plane, which is not acting on space-time, as the only possibility to 
obtain a pure 2-form from a^T* (g) ^ is to have dim{G) = 1. However, we have explained the 
structure of electromagnetism from that of the conformal group of space-time, with a shift by one 
step in the interpretation of the Spencer sequence involved because the "fields" are now sections 
of Ci ~ T* (g) ^ parametrized by Di and thus killed by D2- 

In general relativity, we have similarly proved that the standard way of introducing the Ricci 
tensor was based on a double confusion between the Janet and Spencer sequences described by 
diagrams (1) and (2). In particular we have explained why the intrinsic structure of this tensor 
necessarily depends on the difference existing between the Weyl group and the conformal group 
which is coming from second order jets, relating for the first time on equal footing electromag- 
netism and gravitation to the Spencer ^-cohomology of various symbols. 

Accordingly, paraphrasing W. Shakespeare, we may say: 

" TO ACT OR NOT TO ACT, THAT IS THE QUESTION " 

and hope future will fast give an answer !. 
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